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Abstract 



A macroscopic characterization of fractals showing up a structural transition from 
dense to multibranched growth is made using optical diffraction theory. Such fractals 
are generated via the numerical solution of the 2D Poisson and Biharmonic equations 
and are compared to more 'regular' irreversible clusters such as diffusion limited and 
Laplacian aggregates. The optical diffraction method enables to identify a decrease of 
the fractal dimension above the structural point. 



PACS numbers: 42.20, 05.40. +j, 68.70,+w 



Fractal surfaces observed in nature can become very complex due to structural tran- 
sitions that are generated during their growing processes. Examples of such phenomena 
can be found on bacterial colonies QXJ] and electrochemical deposition experiments Q. 
The observed fractal patterns exhibit a similar class of complex structure despite the fact 
that the mechanisms involved are clearly different. It has been found that, at a certain 
threshold distance, these systems display an intriguing structural transition from dense 
to multibranched growth behaviour of which not much is yet understood. 

In the current literature, two different approaches attempting to explain and repro- 
duce a structural transition during irreversible growth have been proposed. In the first 
one || |], [j| the transition has been derived by solving the Poisson equation on a squared 
lattice -that becomes dependent on the potentials at two boundaries, the distance be- 
tween them, and a screening length. Hence, it has been argued that screening, due to 
free charges, strongly diversifies the patterns that grow in the presence of electrostatic 
fields §. 

An extention of this problem is the second approach due to the authors [0, §], which is 
based on a Biharmonic equation in two-dimensional (2D) isotropic defect-free media. By 
discretizating the Biharmonic equation we have proved that a transition from dense to 
multibranched growth can also be a consequence of a different coupling of displacements 
during the pattern formation. Within the Biharmonic model the transition appears 
when the growth velocity at the surfaces present a minimum as also occurs within the 
Poisson growth. The later implies that both approaches, i.e. Poisson and Biharmonic 
growth, describe a similar class of complex structural transition phenonema from two 
different perspectives and on different systems. 

It is the aim of this work to deal with a macroscopic characterization of fractals 
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showing up a structural transition on growing by using the optical diffraction method. 
The formalism adopted is not new, as it is an extention of what has been largely ap- 
plied to analyse optical diffraction from rather 'regular' sctructures like diffusion limited 
aggregates (DLA) |9| and Laplacian fractals (see, e.g., [[TO], [U], [L2|, |I3[), besides Koch 
clusters [0, [RJ. 

We shall expand on these ideas to study the optical diffraction patterns of fractals 
displaying a transition from dense to multibranched structure. Fractal growth is here 
numerically simulated using the Poisson and Biharmonic approaches and the results are 
compared with the well-known cases of DLA and Laplacian growth. We investigate the 
effect of the structural transition on the diffracted intensity and, consequently, on the 
fractal dimension df for all different types of clusters. 

For completeness we shall give a brief description of the optical diffraction formalism. 
Let us start considering a fractal structure composed of N identical and similarly oriented 
particles on the plane x — y. The position of their centers of mass is given by R n = 
{^niVn)i where n = 1,...,N. In this system we investigate the Fraunhofer diffraction 
pattern for the fractal, assuming that each particle corresponds to one aperture. This 
assumption is experimentally realizable even in the case of dark clusters, since they can 
be photographed on high contrast films which are later lightened Supposing that 
the incident plane wave propagates parallel to the z-axis, the diffraction amplitude due 



to iV 'apertures' can be expressed as |Lq| 



A(k) =C f e- ikr d 2 rJ2e- ikRn , (1) 

J S n=l 

where C is a constant factor, S denotes the surface occupied by one particle. The vector 
k in Eq.(|l]) is the component, parallel to the xy-plane, of the scattered wave vector. Its 
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modulus is given by k = ^sin^ « ^-6, where 9 is the angle (small in practice) which 
the scattered wave vector makes with the z-axis and A is the wave length of the incident 
light. 

The form factor F(k), corresponding to the intensity scattered by one 'aperture', is 
determined by the integral in Eq. (JJ) . The sum over n in Eq. (JJ) determines the structure 
factor 

1 N 

5 ( k )M^Ee- ikR "| 2 • (2) 

iV n=l 

In the case that k is smaller than a -1 (a being the size of an elementary particle), the 
form factor is practically constant, (i.e., F(k) ~ 1) and the light distribution in the 
diffraction pattern is given by the structure factor such that the normalized diffraction 
intensity J(k) m S(k), with ka < 1. 

The factor S(k) of Eq.(|) depends on the distribution of particles in the cluster 
and is related to its fractal dimension. This relation can be found by using density- 
density correlation functions (see, e.g., ||, [13|] and references therein). Accordingly, for 
< k < a^ 1 (L is the size of the whole aggregate) the expected value of the intensity 
-or, alternatively, of the structure factor S(k)- is 

< /(k) >= J d 2 Re" ikR < p(R )p(R + Ro) > , (3) 

where 

< p(R )p(R + R ) >= J d 2 R p(R )p(R + R ) ~ R~ a , (4) 

represents the density- density correlation function. For fractal aggregates this relation 
obeys a power law variation 0, where the exponent a is related to the Hausdorff di- 
mension [17]] df = d — a with d the Euclidean space dimension. This variation of the 



density-density correlation function leads to the power law behavior of the diffracted in- 
tensity as a function of the wave vector I(k) ~ k~ d f 0. Thus, it is possible to estimate 
the fractal dimension df from diffraction patterns. In particular we shall see that this is 
valid for fractal structures displaying a structural transition as those at hand. 

In the following we shall calculate the quantities 5"(k) and < /(k) > of the diffraction 
intensity employing Eqs.(Q) and (|3|) and shall search for a power law behaviour. For 
simplicity we shall assume that the apertures (or particles) have a squared form a x a 
with a distribution obtained by numerically simulating fractal growth. In the simulations 
the particles are distributed on a squared lattice, and the centre-of-mass coordinates x n , 
y n are assumed to take the following discrete values 

x n = m n a, m n = 0, ±1, ±2, y n = t n a , £ n = 0, ±1, ±2, (5) 

Therefore, according to Eqs.(U) and (|5]), the structure factor is calculated by the formula 

1 N 

S(k X ,ky) = | e-™^™^^ | 2 . (6) 

n=l 

Following the DLA simulation model ||, we add particles to the growing clusters 
one at the time undergoing a random walk that starts from a point on a variable circle 
centered on the simulation box. In this process the particles are deposited adjacent 
to occupied lattice sites to then start off again the random walk of a new particle at 
another randomly chosen position, and so on. An example of DLA is shown in Fig.l 
(upper left-corner). 

Laplacian fractals are generated using the growth law (the so-called dielectric break- 
down model) 

V 2 = O . (7) 
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We use lattice sites enclosed within a circle of (normalized) radius r = ^i 2 + j 2 = 100 
such that <p° and <ft l are unity and zero at the outer circular boundary and the inner 
growing aggregate, respectively Seed particles are placed centered in the simulation 
box. The discretazion procedure then follows standard techniques till solutions of 
the discretization of Eq.(|7|) converge to a desidered accuracity. The stochastic growth 
probability P adopted (at the grid site (i, j)) is assumed to be proportional to the local 
field 1, i.e. 
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(8) 



where the sum runs over nearest neighbor sites to a cluster. Herein we set r] = 1 for the 
sake of simplicity. An example of Laplacian fractal is shown in Fig.l (upper right-corner). 

The effect of screening on fractal structures growing under electrostatic fields was 
reported in Ref . || d . Within this model the previuos Laplacian equation is replaced by 
the linearized Poisson-Boltzmann equation 

V 2 <p = \ 2 <f) . (9) 

The origin of screening lies in the presence of free charges and leads to a rich variety 
of patterns -see the example in Fig.l (lower left-corner). The model introduces a new 
length scale, i.e., A, and a nontrivial dependence on the boundary conditions which 
is responsible for a structural transition on growing. The patterns can have a fractal 
character at scales shorter than A, be Eden-like, or grow dense, to then follow the 
transition from dense to single branch growth, which is characterized by a change in the 
sign of the electrostatic field at the aggregate || . However, one of us H has demostrated 
that the transition found using Eq. (|9|) is altered by the existence of a critical field in the 
growth of the pattern. 
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On the other hand, the alternative Biharmonic model under consideration is based 
on the discretization (in a square lattice of size L x L) of the Biharmonic equation JF), H 

V 2 (V 2 m)=0 . (10) 

A crucial difference with respect to the above Laplacian and Poisson models, is that 
within the Biharmonic model, iterative procedures are carried out around thirteen next 
nearest neighbours -and not on four as for Eqs.(^) and (|S|)- on equal grounds. Because of 
this the formation of connected patterns within the Biharmonic equation becomes non 
trivial, hence this model is more involved than Laplacian or Poisson growth. To generate 
Biharmonic fractal patterns -as seen in Fig.l (lower right-corner)- we set for simplicity 
the derivative boundary condition, that is necessary along the radial-direction, equal to 
zero and the growth probability P proportional to V 2 w, (corresponding to the potential 
in Eq.(|)). 

We focus now on the results of diffraction theory. Due to the complexity of the 
present numerical calculations for growing fractals (8 hrs circa of CPU Convex time for 
each run), our results are based in a statistics of only few different clusters for each 
class of complex structures considered. The data shown next is displayed for illustrative 
purposes, whereas the straight lines (used to estimate df) are the result of our crude, 
but representative, statistics. 

To this end we can add that the accuracity achievable depends on the wavevector 
regime considered. Variations of < I(k) > decrease on increasing k and, accordingly, 
the error also decreases. The reason for this essentially lies on the number of particles 
considered (which form the fractal) and also on the shape of the box used in simulations 
(squared in the present case). For k values above the transition point (~ 1/0. 6L) we 
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estimate the error of our calculations to be less than 1%, whereas for k values near to the 
system size ~ 1/L the errors increase to about 5%. Thus, to make a better estimation 
of df it will be necessary to obtain statistic for < I(k) > over a wide range of fc-values. 

As an example we plot in Fig. 2 the diffracted intensity / as a function of the wave 
vector k at the plane k y = 0. The diffraction pattern corresponds to a Biharmonic 
fractal with 1500 particles (each of size a) similar to the one displayed in Fig.l. This 
plot shows strong fluctuations associated to the fractal structure of the aggregate. It can 
be seen that fluctuations between zeros of the form factor (27m <| k x \ a < 2n(n + 1); 
n = 1,2,...) are of the same type. That is, the structure fractor turns out to be a 
periodic function of the wave vector: 

S(k x + 2nm/a) = S(k x ) ; m = 0, ±1, ±2, ... (11) 

This fact is a direct consequence of the assumed distribution of particles on a squared 
lattice (c.f., Eq.(^jD). When the magnitude of the wave vector is sufficiently large, namely 
k x a ^> 1, the intensity decreases considerably due to the form factor F{k x ) ~ k~ 2 . 

Information about the fractal dimension df may be obtained at low values of the wave 
vector ( L^ 1 <C| k x |<C a -1 ). Therein F(k x ) « 1 and I(k x ) « S(k x ). It is noteworthy 
that in the case of deterministic fractals [0 the method, used to obtain the fractal 
dimension df, consists of averaging the structure factor over each of its frequency bands, 
which are scale invariant. The so averaged structure factor < S(k x ) > varies according 
the power law k x d} . In our example the structure factor for a random Biharmonic 
fractal, as shown in Fig. 3, also decreases at k x a <C 1, but it has no scale-invariant 
frequency bands. Therefore, we cannot apply the preceeding relation between S(k x ) and 
k x for random fractals. In order to avoid this problem, we have calculated the averaged 
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intensity denned by the expression 

< I(k) >= — <f d0J(k) ; k x = kcos<p, k y = fcsin0, (12) 

As we shall see next, this quantity has an evident power law behavior as a function of 
the modulus (k) of the wave vector. 

Our numerical results for < I(k) > of diffussion limited and Laplacian aggregates are 
given in Fig. 4 by full and open squares, respectively. The averaged intensity < I(k) > 
is smoother than /(k) and decreases as k~ df . The fractal dimension df is obtained by 
adjusting the log-log plot of < I{k) > to a straight line as indicated in these figures. 
So, the fractal dimensions for DLA and Laplacian patterns obtained from these curves 
are df = 1.72 and df = 1.70, respectively. These df values are closed to the values 
calculated via the box counting method, namely, by counting the number of particles 



N(r) inside an increasing radius r (around a seed particle) -see, e.g., [|l"8[ . 

We shall see next that a structural transition during fractal growth leads to a change 
in their fractal dimension df. The slopes of the curves in Figs. 5 and 6 show the decrease 
of the fractal dimension for Poisson and Biharmonic agregates, respectively. (These were 
roughly estimated similarly to Figs. 3 and 4). 

In the case of Poisson fractals (as the one in Fig.l: lower left-corner), df varies from 
1.60 -full squares in Fig.5- (by counting N=650 particles before the transition) to 1.44 as 
deduced by fitting the open squares to a second line. Before the structural transition the 
df value for a typical Biharmonic fractal (as shown in Fig.l: lower right-corner) is equal 
to 1.66 -full squares in Fig. 6- being quite closed to the corresponding fractal dimension 
of diffussion limited agregates and Laplacian fractals in Fig.4, respectively. After the 
transition point the df obtained from the slope in Fig. 6 (open squares) for < I(k) > of 
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a Biharmonic pattern reduces to a value of 1.59. 

For circular geometry, a very crude analytical analysis of the Biharmonic Eq. (|10|) 
implies that the transition point for Biharmonic fractals is approximately located at a 
distance far from the center equal to r^/(L/2) fa e -1 / 2 by considering displacements of 
the growing surfaces in a continuous limit in reasonable agreement with numerical simu- 
lations |7|, By comparing the structures in Fig.l, it can be seen that Poisson fractals 
display a transition point which appears at a rather shorter distance than irreversible 
Biharmonic structures (the latter enclosed by a circle). This structural transition ap- 
pears at a point in which the growth velocity of the active zone of the clusters exhibits 
a minimum similarly to what occurs within Poisson growth, but with different magni- 
tudes H p!9f] . This fact explains the small differences observed in the value of df after 



the structural transition within the Poisson and Biharmonic approaches. 

Results of optical diffraction thus enables to identify and relate changes in the fractal 
dimension df of aggregates to variations in the diffracted intensity as a function of the 
wave vector. We have seen that there is indeed a decrease of df above the structural 
transition, which appears beyond one half of the system size (assuming t] = 1 in Eq.(|8|)) 
in accord with early crude estimations |7|, ||, ^ made using the box counting method, 
even if it falls beyond the transition point. Different magnitudes of df above the struc- 
tural transitions are due to the different growth velocities obtained from the Poisson and 
Biharmonic approaches. 

Our present findings may -in principle- be experimentally confirmed. The averaged 
intensity < I(k) > might be determined with an experimental arrangement as in the one 
described in |fL0|| . Therein, a photomultiplier is connected to a multichannel analyzer 
to thus record I{k x ,k y ) and the displacement of the photomultiplier is controlled by a 



high-precision motorized micrometer. Then, after scanning the diffraction patterns, the 
average of intensity over concentric circles (< I(k) >) is obtained. 

Concerning our simulations, it would be interesting to extend them to values of 
r\ ^ 1. This is so because by tuning 77 — > 0, the structural transition now corresponds 
to a " dense-to-multibranched transition", whereas for 77 — > 00 one obtains a "transition 
from slow to faster growth" [T9]. 
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Figure captions 

• Fig.l: The different types of fractal aggregates investigated using optical diffrac- 
tion theory. Diffusion limited agregate with 2000 particles. Laplacian fractal 
composed of 2500 particles. Poisson fractal (1300 particles) with structural tran- 
sition. Biharmonic fractal (3000 particles) with structural transition. In the later 
the circle locates the transition point at 60% far from the center. 

• Fig. 2: Central profile (k y = 0) of the diffraction spectrum for a Biharmonic fractal 
with 1500 particles. 

• Fig. 3: Structure factor of a typical Biharmonic fractal as in Fig. 2. 

• Fig. 4: Log-log plot of the angle-averaged intensity < I(k) > against k for diffusion 
limited (full squares) and Laplacian (open squares) agregates with 2000 and 2500 
particles, respectively. 

• Fig. 5: Log-log plot of the angle-averaged < I(k) > against k for a Poisson agregate 
as in Fig.l. Full squares: before the structural transition (650 particles). Open 
squares: after the structural transition. 

• Fig. 6: Log-log plot of the angle-averaged < I(k) > against k for a Biharmonic 
fractal as in Fig.l. Full squares: before the structural transition (1500 particles). 
Open squares: after the structural transition. 
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